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Abstract — Analysis of the effect of the distance between two plates with thick gel layers on the Stern potential 
and disjoining pressure has shown that the interface potential difference inside of the gel layer results in a 
dependence which should be taken into account when further developing the DLVO theory. Analysis of the elec- 
trostatic interaction is based on combining the model of the porous double layer at the oxide—solution interface 
(first introduced by J. Lyklema) and that of charge regulation. Treatment of the electrostatic interaction involv- 
ing the extended model of the porous double layer is of particular importance in the study the stability of bio- 
colloids. Cell glycokalyx is permeable for ions; it has considerable thickness and carries a certain amount of 
ionogenic groups whose degree of dissociation is dependent on the thickness of the electrolyte interlayer and 
on the potential distribution, thus being varied inside of the gel layer. The extended model of the porous double 
layer makes it possible to take into account specific electrosurface properties of the biological cells, contrary to 
the model used by Ohshima, where constancy of the charge is suggested. This fact makes the problem more 


complicated, since potential is varied across the gel layer and calculation of the charge distribution becomes 


necessary. 


INTRODUCTION 


The universally recognized theory of hydrosol sta- 
bility, developed by Derjaguin, Landau, Verwey, and 
Overbeek [1], is closely related to the classical model 
of the electrical double layer for hydrosols. This model 
assumes that the particle surface carries charge, con- 
centrated within an infinitely thin layer, and this charge 
is compensated by the charge distributed over the dif- 
fuse plate. The main parameter characterizing the elec- 
trical double layer in this model is the surface potential 
(¢-potential). When all assumptions of this model are 
fulfilled, the C-potential may be calculated from elec- 
trokinetic measurements, then the value obtained may 
be employed in the DLVO theory formulas to calculate 
the electrostatic repulsion between colloidal particles. 


Comprehensive study of the electrosurface proper- 
ties reveals that the classical model of the electrical 
double layer adequately fits the experimental data for 
most of the actual colloids, and silver iodide sol above 
all. However, similar studies of other colloids, in partic- 
ular, sols of different oxides, demonstrated clearly that 
the classical model fails to correlate measured values of 
different effects. The main discrepancy was found to be 
between very high values of the fixed particle charges 
and the surface conductance of the particles, on one 
hand, and between low values of electrophoretic mobil- 
ity and stability on the other. 


Various corrections to the classical model of the 
electrical double layer have been suggested in order to 
overcome these difficulties. One of the most useful 
hypotheses was proposed by Lyklema in 1968 [2]. 
Lyklema suggested that particles are covered by the gel 


layer, where a significant share of the mobile charge of 
the electrical double layer is concentrated, and only a 
small share of the net charge of the electrical double 
layer (the so-called electrokinetic charge) is distributed 
over the traditional diffuse layer, i.e., over the free solu- 
tion outside the gel layer. This hypothesis was further 
developed by various authors [3 - 6}, being refined by 
various corrections and additional assumptions. Major 
attention, however, was focused by these authors on the 
analysis of the behavior of the surface charge, and 
problems of stability have been only slightly touched 
on within the framework of this model. 


Treatment of the electrostatic interaction in terms of 
the extended model of the porous double layer is of spe- 
cial importance for the analysis of biocolloid stability. 
Cell glycokalix is permeable for ions; it has consider- 
able thickness and carries a certain amount of iono- 
genic groups. It may be assumed that due to these fea- 
tures, the extended model of the electrical double layer 
would result in a better description of the disjoining 
pressure between biological cells than the classical 
model. 


The analysis of the effect of the gel layer on stabil- 
ity, performed by Ohshima [7] for the interaction of 
ionite particles with particles of a different nature 
(including interaction with other ionite particles), was 
the most adequate. Ohshima considered the case of 
constant density of a fixed spatial charge, which is typ- 
ical only for very acidic or very basic ionites. The con- 
stancy of a fixed charge is but a specific case of the gel 
layer models. In general, the charge density is depen- 
dent on the electrical double layer potential. This evi- 
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dence significantly complicates the problem, since it 
involves variation in the potential inside of gel layer, 
thus necessitating the relevant calculation of the charge 
distribution. 


The aim of this paper is to describe a general 
approach to the derivation of the theory of ion-electro- 
static interaction of the particles with charged gel lay- 
ers, and to analyze simple limiting cases of the general 
formulas. 


’ SETTING UP A PROBLEM 


Disjoining pressure in a slit between two identical 
planar surfaces (Fig. 1) is equal to the excess osmotic 
pressure of a solution at the point where solution is not 
affected by the electric force, i.e., in the center of the 
slit. When N types of ions with valence z; and concen- 
tration C; are presented in a solution, then denoting the 
potential in the center of the interlayer with Wo, we can 
express the disjoining pressure P using the following 
formula: 


P = RTY.C(e*-1) (1) 


(from this point on, only dimensionless potential mea- 
sured in RT/F units is used; the symbol X denotes sum- 
mation over types of ions from 1 to N; index i indicates 
a definite type of ions). 

The potential profile in a double-layered medium 
consisting of both the solution interlayer and a suffi- 
ciently thick gel layer should be derived in order to cal- 
culate Wo. Consider first the boundary conditions of this 
problem. By denoting the permittivities of the solution 
and the gel by €, and €,, respectively, we achieve the 


values 
4nF” ‘i 
Kae = | RTE, oe @) 


which are inverse screening lengths in these media. 
Coordinate x is normalized to the screening length in 
the respective medium: 


X,= KX; X, = KX. 


Making use of these notations, the boundary conditions 
may be written as 
dw 


es = 0, (3) 
dx,| , 


,=0 
which is a condition of the symmetry of the potential 
distribution; 

Viz,2n-0 = Wis =nso (4) 


is a condition of the potential continuity at the phase 
boundary (if X, is the total interlayer width, then 
No. 6 
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Fig. 1. Coordinate system and notations to figures. 


H = x,X, /2 is the dimensionless half width of the inter- 
layer); 

d 
= KE . 


as 
x,=h-0 dx, 


ay 
Ke say, 


(5) 


x, =h+0 
acondition of the continuity of the dielectric induction; 


dy 
mote = 0, 6 
dx, XxX, 0 
and 

Mo sa = Mp (7) 


are conditions of the uniformity of the gel volume; 
Wp is the potential of the bulk phase of the gel. In the 
theory of ionites, this value is known as the Donnan 
potential. This potential is independent of any inter- 
layer properties and is solely determined by the proper- 
ties of both the gel, and the solution which is in equilib- 
rium with the gel. Since the uniform volume is neutral 
electrically, the potential yp may be determined from 
the equation 


P.(Wp) = 0. (8) 
A specific expression for gel charge density p, will be 
discussed further. 


In solution, the potential obeys the Poisson—Boltz- 
mann equation 


dy ~z; Cz; 
a “yes 4, = Foe o 


In gel, the ion concentration is not expressed in 
terms of potential through the Boltzmann distribution; 
first, the density of active sites and, hence, the maxi- 
mum ion concentration, are finite; second, ions interact 
with active sites in a specific way, i.e., not electrostati- 
cally. Two variants of distribution are considered most 
often. The first case represents independent adsorp- 
tion [4] when there are adsorption potentials ®, for each 
type of ion, and site densities Q; are attributed strictly 
to ions of the ith type. The second case refers to com- 
petitive adsorption [5S], when net site density is intro- 
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duced and each single site may be occupied by an ion 
of any type or remain free. Theoretically, we may also 
consider any intermediate case which may be a combi- 
nation of these two variants. A number of problem 
parameters, related to the ionic components, may range 
from N, for competitive adsorption, to 2 x N, for inde- 
pendent adsorption, respectively. Here, we focus atten- 
tion on the independent adsorption model. Concentra- 
tions of the adsorbed ions A, are 


Ps gad 
A, = C,B;; “B; = ——.——_, (10) 
1+ Ci O26 
Q; 
and the charge density is 
P, = FY zA,. (11) 


Expression (10) is, in fact, an adsorption isotherm of 
the Langmuir—Stern type; this distribution is converted 
to the Boltzmann distribution at ®, = 0 (no specific 
interaction between ions and gel) and Q; —> © (no sat- 
uration). Potential within the gel layer obeys the Pois- 
son equation where charge density is expressed by for- 
mula (11): 


(12) 


In compliance with [2], this equation will be 
referred to as the Poisson—Langmuir equation. 


GENERAL SOLUTION 


Since the right-hand terms of equations (9) and (12) 
are dependent solely on the unknown function, one may 
readily obtain the general solution of these equations in 
quadratures. Single integration of the equation (9) with 
account for (3) yields the following solution: 


dy \ r,; Z 
—'| =2) —e =e"): (13) 
( dx, py qj 
Double integration of equation (13) gives 
V, ro 
h= (14) 
I 


Yo {2 “ig ay Vo 
> Z (e e'’) 
For the potential yp, we transform equation (8) to the 
following form: 


—O,-7Vp 


ee se 
{eater 


i 


(15) 
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The single integration of equation (12), with account 
for equations (6) and (7), yields 


C; ~O,-7M, 
Vv 2 err} 
=i = 2))v;In moe 
9, 


pera 


n 1a ~O,- 2M, 
i 2M, —ZVo: ies i 
i l+—e ' 
Q; 
(17) 


Having calculated yp by solving equation (15), we 
substitute the resultant expression into equation (17). 
Then equations (14) and (17) form a system of two 
transcendental equations needed to calculate y, and Wo. 


Maximal disjoining pressure in our case is achieved 
at surface contact, i.e. at h —» 0. As seen from (14), in 
this case Wo —> Y,, i.e., the potential inside a very thin 
layer remains constant. Substituting Wo = y,, into equa- 
tion (17), we obtain the solution of this equation: 


Wh = Wo- (18) 
Furthermore, inserting equation (8) into (1), we get 


Prox = RTD Che“? - 1). (19) 


This expression is valid for the considered case of infi- 
nite gel layers, as well as for the interaction of gel lay- 
ers of arbitrary thickness. Indeed, since at the coales- 
cence of two gel layers a spatially uniform phase is 
formed, the potential of this phase is bound to take the 
value Wp, which is independent of geometrical param- 
eters of the phase, such as its thickness. 


The fundamental importance of formula (19) is that 
it indicates the quantitative difference between the 
behavior of the particles, with and without the gel layer. 
The electrostatic component of the disjoining pressure 
for the particles without gel layers is solely governed by 
the surface potential, whereas for the particles with gel 
layers the controlling factor is potential wp, which is 
always higher than the surface potential. Since any 
complete description of the particles with gel layers 
requires knowledge of potential yp, a measurement of 
only ¢-potential cannot provide all the necessary infor- 
mation. 


LOW POTENTIAL DIFFERENCE 


Hence, the potential difference between the depth of 
the gel layer and the middle of the infinitesimal thick- 
ness interlayer is lacking. Therefore, at small thick- 
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nesses this difference is small. We take advantage of 
this fact to obtain an approximate solution. We denote 


AW, =Wo- V3 AW2=W,-— Wo; [AY] <1. 


Expanding integrand of expression (4) in the vicin- 
ity of the lower limit, we obtain 


|AV, 2M 
h= 8B,’ Bo = -Sir-e ‘ 


Within the framework of this approximation, equa- 
tion (17) acquires the form 


€,ByAy, = €,B,AW;, 


(20) 


(21) 


[When expanding right-hand part of equation (17), the 
term linearly dependent on Ay, vanishes, in view of 
equation (15).] 

Substituting equation (20) into (21), we obtain 


Bee oN 8Byh 
E, B, ? We 7 0 ° 
Hence, at h —> 0 inequality |Ay,| > |Aw,| is valid and 
virtually the entire potential difference is concentrated 
inside of the gel; ignoring the terms of the order 
of O(h?), we finally get 


; 8E, 
Wo = Wp- Bo €B, 


Calculating By from formula (20), we may substi- 
tute Wp for Wp. This substitution results in an error of 
the order of O(h”) in equation (23). The only limitation, 
while using equation (23), is that Wo should not be too 
low since at [Wl = 0 Bo vanishes, and formula (20) 
should be modified. In this case, the Debye approxima- 
tion should be employed. 


Ay, = Boh (22) 


(23) 


‘THE DEBYE APPROXIMATION 


Potential in the system studied is monotonically var- 
ied from Wo to Wp. Therefore, if 


IWpl <1, (24) 


then the potential at any point of a system will be low 
in absolute magnitude. Using this fact, one may readily 
determine integral (14) as 


Wi), = Wocoshh (25) 
and calculate the Donnan potential using equation (15): 
C; -« 
Wp = B./B;; “9° 5 
VEZ vieizi o 
Lise? Bho yey 
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Expanding equation (17) into a series by the powers of 
entering potentials [with accuracy up to O(h*)], and 
making use of equations (25) and (26), we obtain the 
following quadratic equation for y,, 


BY, = Vip — 1) - 2,4 + ¥58 — Wp), 
2¢, tanh*h 

ae ~€,B; 

with only the smaller value of the square root having 


physical meaning 
1+B, ° 


because only this root fits the condition y, —> Wp 
ath —> 0. 


(27) 


bY 


(28) 


W, = —-Wo 


APPLICABILITY LIMITS 
OF THE APPROXIMATE FORMULAS 


Problems of the applicability limits for the two 
approaches discussed above have been specially treated 
in [8]. It was established that for 1:1 electrolyte it may 
be assumed (with 10% relative error) that an approxima- 
tion of a small potential difference is valid at h < 0.45, 
with the Debye approximation at |Wp| < 0.57. We 
believe that these estimations may be taken as reference 
points, even in more complicated cases than those 
treated in [8]. 


DEGENERATION 
Analyzing the expression (10) for the concentration 
of the adsorbed ions, we may note that this expression 
is readily degenerated, i.e., it becomes more simple in 
the sense that its dependence on some input parameters 
ceases to exist. Then, the key value is 


C; -O -zW 
_=—e i ‘i . 
uO, 


It may be recognized that, in a nontrivial way, expres- 
sion (10) is affected by all entering parameters only at 


0.1 <x; < 10. (30) 


(29) 


At 

xX; <0.1, (31) 
coverage of the active sites is so small that the role of 
coverage may be neglected (the Henry limit) and 


A; = Cee, (32) 


As can be seen, in this limiting case the dependence on 
Q, disappears. In the opposite case, at 
12 10, (33) 


coverage of the sites is close to unity (the saturation 
limit) and, hence, coverage becomes constant. Then 


A;= Q;, (34) 
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that is, dependence on C,, ®,, z;,, and w ceases to exist. 

The narrower the experimental range of parameters, 
the higher the probability that a part of the ions will be 
susceptible to degeneration. In this case, the restoration 
of a set of parameters from experimental data is 
strongly complicated due to the fact that we are trying 
to evaluate a larger number of parameters than are actu- 
ally involved in the problem. Degeneration leads to 
increasing dispersion of the evaluation of the mostly 
degenerated parameter, and of all other parameters cor- 
relating with that parameter, because a larger number of 
parameters is responsible for the same bulk of informa- 
tion. Multivariate algorithms, including those involved 
in the analysis of the significance of all estimated 
parameters, should be employed in order to obtain the 
most reliable estimations. 


SIMPLE PARTICULAR CASE 


All of the calculations performed are valid in the case 
of an arbitrary N-component solution which is character- 
ized by 2 + 2N parameters. To get better understanding 
of some features of our problem, consider the particular 
case of a 1:1 binary electrolyte. In this case, 
ral =+1 »2=-l ’ C,=C,=C, x, = +0.5, A, =-0.5. 
For further simplification, assume that ion 2 is inert 
(® = 0) and ion 1 is potential-determining (®, < 0, 
|@,| > 1). Then equation (15) for the Donnan potential 
will be converted to 

~®,-Vp Vp 
= Cen Tow ee 
jee ee 


Q, Q, 
with 
Wp =-®,/2+ In(Ji+q° +4), 
- C(Q, - Q») oo? 
22,2, 


Even this example [equation (36)] demonstrates the 
degeneration of the problem: i.e., the dependence of yp 


Wo 
20 


(36) 


10E0 


10E-2 
10E-5 C, mol/l 


10E-6 


Fig. 2. Dependence of the Donnan potential of the gel layer 
on the electrolyte concentration (Q; = 10): ( - 3) ®, =-1, 
(4-6) ®, = -10; (1, 4) Q, = Q,/10, (2, 5) Q = Q,, 
(3, 6) Q2 = 10Q;. 
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on a combination of Q, and Q, rather than on Q, 
and Q,, separately. 

At Q, = Q, (the density of the adsorption sites is the 
same for both types of ions), 


Wp =-®, /2 (37) 


at any electrolyte concentration. The same situation 
will be observed at any values of Q, and Q,, when the 
electrolyte concentration is sufficiently low (the Henry 
limit, g < 1). Under saturation at sufficiently high elec- 
trolyte concentration, 


Wo =—®, + In(C/Q,); Q, < Qy, (38) 
Wo=1In(Q,/C); Q.> Q). (39) 

In either case, the degeneration “takes up” one or two 
parameters. When Q, # Q,, it is especially difficult to 


determine the larger value against the background of 
the smaller. 


When Q, = Q2, formula (23) is simplified to 


x 8e 
W p= —®,/2 + hsinh(®,/2)e" "(1 + se a4) —. 
1 
(40) 


Figure 2 shows a few plots calculated by equa- 
tion (36), illustrating the dependence of the Donnan 
potential on concentration. At Q, # Q3, each curve plot- 
ted in semilogarithmic coordinates may be distinctly 
divided into three, regions: two straight lines corre- 
sponding to limiting cases (38) and (39), connected 
together by a narrow (smaller than one decimal order of 
magnitude) transition zone. When the adsorption 
potentials are low (®, = —1, in Fig. 2), the position of 
the transition zone is determined as C = min(Q,, Q,). 
In Fig. 2, the degeneration is clearly visible, both at low 
and high concentrations [at Q, = 10Q, and C> 10“ mol/l, 
the curves for yp coincide irrespective of @,, due to 
condition (39)]. 

Figures 3a and 3b demonstrate data calculated by 
formula (40) for €, = 15. At very low concentrations, 
the disjoining pressure is small due to the very low con- 
tent of ions in the system [see formula (1)]. When the 
concentration is raised, pressure first increases 
(the ascending part of the curve in Figs. 3a and 3b), 
and then starts to reduce rapidly due to a decrease in 
the overlapping of the electrical double layers, and a 
sharp decrease in Wo. This is qualitatively consistent 
with the pattern of the disjoining pressure for the 
model of a smooth surface with constant surface 
potential. 

Differences in these models can be followed with 
the help of the Debye approximation (Fig. 4, Q, = Q,, 
€, = 15). In this figure both surfaces are characterized 
by a constant potential equal to the Donnan potential of 
the gel. As a result, both models appear to be alike, pro- 
vided that C < Q,. In the case of a smooth surface, the 
potential in the middle of the interlayer reduces solely 
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Pd (a) Px10°,Pa 


5.0E-7 C-inol/l 


Fig. 3. Dependence of the disjoining pressure between two 
gel layers on the electrolyte concentration, calculated in the 
approximation of low potential difference using equation we 
Q; = Q) = 10% mol/l; clearance width (X4) = 20 A; 
(1) dimensionless pressure (on the left-hand axis Y), (2) pres- 
sure, in Pa (on the right-hand axis Y); (a) ®; =—1, (b) B, = 
—10, (1) dimensionless pressure, Pa, (2) pressure P, in atm. 


due to the decrease in thickness of the electrical double 
layer, whereas in the case of gel layers, rearrangement 
of the potential between gel and solution, accelerating 
the reduction of the potential with increasing concen- 
tration, is also important (in addition to this first factor). 
This circumstance is revealed by the fact that the steep- 
ness of the descending part of the curve for dimension- 
less disjoining pressure is higher for the gel layer and 
lower for the smooth surface. When capacity Q, is 
increased, curve / for the gel layer shifts to the top right 
of the curve until it coincides with curve 3. With 
increasing interlayer thickness, curve 3 shifts to the left, 
whereas the slope of the curve J (in addition to the shift 
to the left) decreases until both curves will merge again 
into one curve. At low potentials, the disjoining pres- 
sure of gel layers will always be lower than pressure of 
smooth surfaces, with the C-potential equal to the Don- 
nan potential of gel layers. Conversely, the disjoining 
pressure of gel layers will always be higher than that for 
smooth surfaces, with the C-potential equal to C-poten- 
tial of the isolated gel layer. 


CONCLUSION 


The main parameter controlling the electrostatic 
component of the disjoining pressure of the particles 
covered with the gel layer is the Donnan potential of the 
gel. This value is beyond the reach of direct measure- 
ment. However, it may be calculated by complex elec- 
trosurface measurements, including the detection of the 
surface charge by titration and the measurement of sur- 
No. 6 
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Pd Px 10°, Pa 
0.15 “41B-1 
] 1E-2 
1E-3 
1E-4 
1E-5 
1E~6 
J 1E-7 


0.10 


0.05 


1.0E-6 
1.0E—-5 
C, mol/1 


Fig. 4. Dependence of the disjoining pressure between two 
gel layers on electrolyte concentration, calculated in the 
Debye approximation by equation (36): (J, 3) dimension- 
less pressure (on the left-hand axis Y), (2, 4) pressure, in Pa 
(on the right-hand axis Y);-(/, 2) gel layer model, 
(3, 4) smooth surface model; (/) dimensionless pressure P, 
(2) pressure P, in atm. Q; = Q2= 10 mol/1; clearance width 
(X4) = 20 A; ®; =-1. 


face electrical conductivity. When the Donnan potential 
is replaced by the C-potential, which can be measured 
directly ysing electrophoresis, significant error may 
arise during the calculation of the disjoining pressure. 
This circumstance is very important when studying 
biological cells, where molecular attractive forces are 
low and, hence, where small variations in electrostatic 
repulsive forces may dramatically change the total 
force balance. 


Three approximate methods, employed to receive 
analytical results, may be applied both to the case of 
independent adsorption (which was treated in this 
paper) and to that of competitive adsorption. These 
include expansion over the small potential difference, 
the Debye approximation, and the approximation of the 
thin double layer (see Appendix). 


Any variant of the description of the fixed charge of 
colloidal particles by means of the Stern—Langmuir iso- 
therms involves a large number of parameters. Theoret- 
ical formulas used to estimate these parameters from 
experimental data are susceptible to degeneration, and 
thus, parts of these parameters defy possible estima- 
tion. This difficulty may be overcome by using ideas of 
experimental design. 
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APPENDIX 
Case of a Larger Width of Interlayer 


When h —» 9, the potential in the interlayer center 
tends to zero. We denote 


We = Wilaoew? AY, = We Ve 


where y; is the potential of the gel—solution interface 
for a single gel layer. w,Potential fits the equation 
resulting from (17): 


Aes Se 
rag as 1) 


(Al) 


(A2) 


Expanding (17) over powers of small values of Wo 
and Ay,, and with regard for (A2) we obtain 


C; —®,- 74, 


€,(AW3 — Wo)/2 = AW3€,)/——— 
l+—e ' 
Q; 

As can be seen from this equation, Ay; ~ OCW): Lé:, 
at very slight overlapping of the electrical double layers 
in the interlayer, the potential of the gel surface varies 
negligibly in comparison to that in the interlayer center. 
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Hence, with good accuracy one may calculate wy from 
equation (14), substituting y, for w;: 


vy 


dy 


| ea re 
J hes ase Wo 
: 2 SMe -y-%| 


Here, the simplification is that the system of equa- 
tions (14) and (17) is decomposed into two indepen- 
dent equations: equation (A2) to determine wy and 
equation (A3) to determine Wy. 


(A3) 


According to [8], this approximation is valid at an 
error level of 10%. 
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